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1.  Introduction 


This  paper  addresses  the  following  inverse  problem:  Can  one  determine 
an  unknown  conductivity  Y  inside  a  body  Cl  by  means  of  static  measurements 
at  the  boundary?  Since  there  is  no  dependence  on  time,  the  underlying  equation 
7*(yVu)  *  0  is  elliptic.  Our  first  paper  [6]  established  the  identifiabilitv 
of  a  real-analytic  y  ;  in  this  one  we  demonstrate  a  similar  result  when  y 
is  piecewise  real-analytic  (for  example,  piecewise  constant). 

Consider  a  bounded  domain  Q  €  Hn  ,  n  z  2  ,  and  a  "conductivity" 

oo  J 

Y  €  L  (0)  ,  0  <  Yq  S  y(x)  .  For  appropriate  $  there  is  a  unique  u  €  H  (Cl) 
such  that 


(1) 


7*(y(x)7u)  *  0  in  SI 

u  “  on  3£2  , 


obtained  for  example  by  minimizing  Dirichlet's  integral.  Let  Q ^(t>)  be  the 
energy  of  the  solution. 


qy<*)  - 


If  3 H  is  sufficiently  smooth  then  Green's  formula  gives 


Q  (*) 
Y 


J  3  a 


3u 

u’Y  —  ds 

3v 


where  "ds"  denotes  surface  area;  consequently  Q  (f>)  depends  only  on  the 


Dirichlet-data  $  and  its  associated  Neumann-data ,  the  conormal  derivative 
3u 


Y  —  .  Our  inverse  problem  is  t£  determine  y  given  knowledge  of  the  quadratic 
.  We  shall  say  that  y  is  identifiable  (within  a  certain  class) 


3  v 

form  Q 


bv  boundary  measurements  if  the  map  y  "*  Q  is  injective  (in  this  class). 


There  is  an  equivalent  formulation  in  terms  of  the  map  taking  Dinchlet- 


to  Neumann-data 


A  •  U  -  y  — 

Y  ‘3ft  Y  3v‘ 9a  ' 


The  form  determines  the  map  ,  by  polarization.  Therefore  our  problem 

is  alternately  to  determine  y  given  knowledge  of  the  Dirichlet-  to  Neumann-data 


map  A^  . 

3y 

If  y  and  •—  are  known  at  3ft  ,  then  there  is  another  formulation 
1  3v 

1/2 

involving  -  A  +  q  .  Indeed,  v  *  y  u  solves 


-  Av  +  qv  *  0 


in  ft 


1/2. 

v  *  y  on  3ft  , 


1  IVyI  ,  1  Ay 
"4  y2  *  2  Y 


whenever  u  solves  (1),  assuming  sufficient  regularity  of  y  .  The  Dirichlet- 


to  Neumann-data  map  for  q  is 


r  I  3v , 

Aq  :  VUft  3  v  1  3  ft  * 


If  the  boundary  values  of  y  and  —  are  known  then  information  about  A 

3v  q 

is  equivalent  to  that  about  A  ;  hence  our  problem  is  also  find  the 

unknown  potential  q(x)  given  knowledge  of  the  map  A^  . 

Some  special  cases  of  the  inverse  problem  can  be  treated  by  separation 

of  variables.  R.  E.  Langer  considered  a  iavered  half-space  in  1933  [8],  He 

reconstructed  all  the  derivatives  of  y  at  the  boundary  from  knowledee  of  A  V.'  ) 

>  J 

for  just  one  particular  flux  •  Moreover  his  choice  of  , ^  is  one  which  natural lv 


shows  up  in  geophysical  applications.  Another  example  is  that  of  a  layered 
rectangle,  analyzed  by  Cannon,  Douglas  and  Jones  in  1963  [2]. 


The  general  case,  without  restrictions  on  the  form  of  ft  and  y  ,  has 
only  been  considered  more  recently;  to  the  best  of  our  knowlege  it  was  first 
raised  by  A.  P.  Calderon  [1].  He  proved  that  the  map  y  is  Frechet 

GO 

differentiable  for  y  Z  L  ,  and  that  th$  differential  at  y  =  constant  is 
injective.  However  its  range  is  not  closed,  so  the  implicit  function  theorem 
does  not  apply,  and  one  cannot  conclude  the  identif iability  of  y  from 
this  analysis. 

We  took  a  different  approach  in  [6],  using  the  variational  principle 
and  special  Dirichlet  data  with  localized,  highly  oscillatory  behavior.  We 
proved  that  the  quadratic  form  (or  the  map  determines  all  the 

derivatives  of  y  at  the  boundary.  Identif iability  in  the  class  of  real- 
analytic  y  follows  as  an  immediate  corollary.  A  more  extensive  review  of 
the  literature  is  found  in  [7],  It  seems  worth  noting  that  the  method  of  [6] 


can  also  be  used  for  -  A  +  q  ,  even  when  q  does  not  have  the  form  (2).  If 

00  ** 
q  and  3ft  are  C  then  the  Dirichlet-to  Neumann- data  map  or  the 

2  2 

corresponding  energy  form  q  /(|vv|  +q|v[  )dx  (provided  each  is  well-defined) 


determines  q  and  all  its  derivatives  at  the  boundary. 


Unfortunately,  the  results  just  summarized  are  far  from  satisfactory. 


The  original  goal  was  to  reconstruct  y  everywhere,  while  these  results 


give  only  its  derivatives  at  the  boundary.  It  remains  unknown,  for  example. 


whether  a  C°°  conductivity  is  determined  in  the  interior  by  these  boundary 
measurements. 

The  present  paper  represents  some  modest  progress  toward  the  interior 
reconstruction  problem.  Our  main  result  is  one  of  identif iability  in  the  class 
of  piecewise  real-analytic  conductivities.  The  proof,  presented  only  in  dimension 
two  for  ease  of  exposition,  uses  the  ideas  of  [6]  together  with  the  Runge 


Approximation  Property  for  solutions  of  V*(y7u)  «  0  .  Its  concept  is 

sketched  in  section  2,  and  the  details  are  developed  in  sections  3-5.  Our 

second  result,  presented  in  section  6,  is  the  identif iability  of  a  layered 

conductivity  which  is  merely  three  times  continuously  differentiable .  For 

2 

simplicity  we  treat  only  the  case  of  a  finite-width  strip  in  ]R  .  The  proof 
uses  separation  of  variables  and  recent  progress  in  one-dimensional  inverse 
spectral  theory.  Our  third  topic,  presented  in  section  7,  is  the  convergence 
of  £  reconstruction  algorithm  for  real-analytic  y  .  This  algorithm  restricts  its 
attention  to  approximations  of  y  within  a  suitable  compact  set  and  the 
proof  of  convergence  depends  heavily  on  the  identif iability  demonstrated  in 
[6].  The  proof  gives  no  information  on  the  rate  of  convergence  or  the 
efficiency  of  the  algorithm; these  remain  important  directions  for  further  study. 


We  begin  by  reviewing  how  boundary  measurements  were  shown  to  determine  a 
real-analytic  conductivity  in  [6].  Given  a  smoothly  bounded  domain  7.  ,  consider 
two  coefficients  and  ,  and  let  u^  and  u?  be  the  corresponding  solution 

of  (1): 

(3)  V*(viVui)  =  0  ,  ui | =  *  ,  i  =  1,2  . 

If  then  their  Taylor  expansions  differ;  relabeling  if  necessary, 

one  easily  concludes  that 

Y1(-)"Y2(i)  '  .  P(x)  =  dist  (x, 3fi) 

for  some  t  >  0  ,  in  a  71  -  neighborhood  D  of  some  x^  €  3f2  .  By  choosing 
Dirichlet  data  $  that  oscillate  rapidly  near  x^  ,  one  can  arrange  that  the 
energy  of  u^  be  concentrated  in  D  ,  and  indeed  that 

(4)  I  |Vu  j2dx  <  z  [  p'*'|Vu1|2dx 

h\D  1  '  D 

with  e  >  0  as  small  as  desired  (see  Lemma  3) .  Then 

2  2 

(5)  Q  (<t>)  =  Y,|Vu  |  dx  >,  Y,  |Vu -|  dx 

Y1  Jq  1  1  JD 

2  Z  ? 

>  Y2  I  i  dx  +  C  P  j  7u^  |  "dx 
D  D 


>  I  Y,|Vu  |2dx 
hi  ~ 

’  2 

(6)  >  Yt  |?U  I  dx  =  Q,  (d>)  , 

a  ‘  r2 


using  (4)  and  the  variational  principle  for  (3)  in  the  last  two  steps.  Therefore 


for  this  choice  of  <t>  the  boundary  measurements  are  different. 


Our  method  for  the  piecewise  analytic  case  is  a  direct  extension  of 
this.  Consider  the  simple  example  of  the  ball  3  =  {|xjcl}  with  an  unknown 
concentric  inclusion: 


Y.  (x) 


lii!<ri 

1*1 >ri 


i  =  l,; 


Here  y|  and  yV  are  real-analytic  functions  of  x  .  If  y^  and  Y 2 
b  a  the  same  boundary  measurements  then  the  argument  sketched  above  shows 
that  y”  =  Y^'  in  i 2S !  >  maxCr^,^)  .  A  new  idea  is  required,  however,  to 
establish  that  y^  =  y^  also  for  |xj  <  max(r^,r9)  . 

Assuming  that  ri  '  r2  ’  w^at  we  wou^^  really  like  to  do  is  to  repeat 
the  estimation 


1 2  '  1 2 

y. I Vu  |  dx  >  I  Y? | Vu  |“dx 

JQ  JQ 

with  n  replaced  by  the  subdomain  3'  =  { |x[ <r? }  .  At  first  glance  this 
seems  impossible,  since  the  elliptic  equation  V»(yVu)  =  0  will  smooth  out  any 
oscillations  of  the  Dirichlet  data  prescribed  on  the  outer  boundary.  But  in 
fact  it  is  possible,  since  this  equation  has  the  Runge  Approximation  Property: 
if  V»(yVu)  =0  on  a  subdomain  u  c  C  ,  then  u  can  be  approximated  on 
compact  subsets  of  00  by  solutions  of  the  same  equation  in  the  full  domain 
(see  Lemma  2). 

Therefore  to  show  that  y^  =  y,,  across  jxj  =  r9  ,  we  once  again  assume 


the  contrary:  if  not,  then 


for  some  (  ^  0 


,  in  a  neighborhood  D1  relative  to  S'  =  {jxj<r0i  of 
some  Xy  *  38'  .  We  suppose  for  simplicity  that  y^  is  the  larger  in  this  neighbor¬ 
hood.  The  construction  oi  [6]  together  with  the  Runge  property  leads  to  a  function 
u^  with  v*(y^7u^)  =0  on  all  of  ~  and 

jlu  j“dx  <  z  [  (c*)^|7u  j^djc 

Jn'\D'  1  -fD'  1 

(see  Lemma  4).  The  behaviour  of  u^  off  JJ'  is  essentially  unknown,  but 

that  is  acceptable  because  y^  =  y^  there.  Estimates  parallel  to  (5 ) — ( 6 ) 

using  <(>  =  u  |  show  that  Q  (4>)  >  Q  ($)  when  e  is  small  enough. 

1  Y,  Y 2 

Therefore  if  y^  and  y^  had  the  same  boundary  data  they  must  have  been 
equal  in  a  neighborhood  of  !  =  r?  . 

If  r^  =  r 2  this  argument  establishes  y^  =  •  If  r^  <  r 2  then 

we've  shown  that  y^  =  y ^  in  {|}t|>r^}  ;  the  same  argument  can  be  repeated 

at  \x\  *  r^  to  establish  y|  =  y^  in  {  |  sc  |  <r^ }  ,  and  it  follows  that 
y^  =  y^  .  Though  the  geometrical  and  analytical  technicalities  are  substantial, 
this  procedure  of  "marching  inward  from  the  boundary"  will  work  in  essentially 
the  same  way  for  the  general  piecewise  analytic  case. 
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3.  Preliminaries 

We  begin  with  definitions  and  conventions  on  the  use  of  the  term 

"piecewise  analytic".  From  this  point  on  we  drop  the  prefix  real  and  use  the  term 

ana] vt ic  as  a  synonym  for  real-analvt ic .  Let  w  he  an  on on  set  in  TR™  . 

m  =  1  or  2  .  We  recall  that  a  function  y  :  us  -*■  TR  is  analvt ic  in  jj  if  it 

is  infinitely  often  differentiable  and  for  every  compact  subset  K  c  there 

exists  a  constant  r  for  which 

K 

sup  rK°^  <  “  * 

Che  supremum  being  over  all  points  x  €  K  and  all  m-tuples  a  of  nonnegative 
integers. 

We  shall  say  that  a  C°°  function  y  is  analytic  on  ui  if  it  has 

an  extension  which  is  analytic  in  a  neighborhood  of  u>  .  A  mapping 

2 

g  :  [0,1]  -*■  E  is  called  an  analytic  curve  if  each  of  its  coordinates  is 

analytic  on  [0,1]  and  Dg(s)  ^  0  for  all  s  €  [0,1], 

2 

A  bounded,  open  set  u>  c  is  called  a  piecewise  analytic  domain  if 

(i)  to  is  connected, 

(ii)  its  boundary  9u>  is  a  union  of 
finitely  many  (images  of)  analytic 
curves , 

(iii)  oi  lies  locally  on  one  side  of  . 

Notice  that  the  boundary  of  a  piecewise  analytic  domain  can  have  cusps. 

*> 

Following  [3]  we  say  that  a  bounded,  open  set  j  c  R“  is  an  analytic 

curvilinear  polygon  if  it  is  connected,  and  if  for  each  x  ]  there  exists 

2 

x  (relative  to  ]R  )  and  a  map  $ 


a  neighborhood  B  of 


such  that 


(i)  $  maps  B  injectively  onto  a  neighborhood 

2 

of  the  origin  in  F*" 

(ii)  $  and  $  ^  have  analytic  coordinate  functions 

(iii)  4> (ojDB)  is  either  tx  :  x2>0;  0  -t>(B)  ,  \x  :  >  0  and 

xo^0-  P'5)  or  -x  :  x  -0  or  x9>0]  0  $(B)  . 

Analytic  curvilinear  polygons  are  special  cases  of  piecewise  analytic  domains 

without  cusps  and  "corners’1  of  angle  t  . 

N 

A  family  {oi.K  ^  is  a  (disjoint)  piecewise  analytic  cover  of  a  closed 
set  co  if 


(i) 

each 

01 

J 

N 

(ii) 

Z  £ 

j 

1 1 

=  1 

(iii) 

o>^  n 

01  . 
J 

Y  :  w 

-  R 

is 

A  function  y  :  ui  R  is  piecewise  analytic  on  oi  (relative  to  the  cover 

{ u> . } )  if  it  is  analvtic  on  each  oi.  ,  1  $  j  <  N  .  Notice  that  no  continuity 

J  j 

is  assumed  across  the  interfaces. 

In  proving  identif iability  we  shall  consider  two  conductivities  y^  anc 
y2  corresponding  to  different  piecewise  analytic  covers.  It  is  important  to 
note  that  they  are  both  piecewise  analytic  relative  to  a  common  refinement  of 
the  two  covers; 


Lemma  1.  If  y ,  and  y  are  piecewise  analvtic  functions  on  oi  ,  then 


there  exists  a  oiecewise  analvtic  cover  {oi..l\  ,  such  that  both  y,  and  y 


are  analytic  on  each  oi^  ,  1  £  j  £  N  . 

(i)  Ni 

Proof:  Let  { uj v  '}  ,  be  a  cover  relative  to  which  y.  is  piecewise  analvtic 

-  j  j=l  l 

(1)  -  (2) 

i  =  1,2.  For  any  pair  (j,k)  the  boundary  of  the  intersection  ^  oj, 
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consists  of  finitely  many  (images  of)  analytic  curves.  Of  course  fl 

1  k 

may  not  be  connected;  also,  as  shown  in  figure  1,  the  boundary  might  intersect 


a)  „  (2) 
JJ  rH' 


Cl \  -> 

=  {jcfR” :  1/2<:_xI<1j 


=  {x€E"  :  lx- (0,1/4)  1  c  3/4  1 


Fig.  1 

itself  at  isolated  points.  Nevertheless,  by  considering  the  connected  components 
and  "cutting  out"  a  neighborhood  of  any  point  of  self-intersection  and  dividing 
this  neighborhood  in  separate  pieces,  wo  can  obtain 

'  (1)  o  (2j  ”  ~ 

W  .  I  I  W  -  J  oJ  , 

J  R  «-l  1 

where  the  u*  are  disjoint,  piecewise  analytic  subdomains  of  ^  • 

The  collection  of  all  w*  (corresponding  to  all  pairs  (j,k))  constitutes 

v 

a  cover  of  Q  relative  to  which  both  and  y?  are  piecewise  analytic. 


Remark  _1 . 

In  the  preceding  proof  we  used  the  fact  that  two  analytic  curves  either 

intersect  in  at  most  finitely  many  points  or  the  intersection  is  itself  an  analytic 

curve.  This  was  used  to  conclude  that  the  intersection  of  two  piecewise  analytic 

domains  has  a  boundary  which  consists  df  finitelv  many  analytic  curves.  A 

k. 

similar  fact  does  not  hold  for  C  -curves,  and  this  is  why  we  work  with  piecewise 
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analytic  part  ions  of  ft  .  . _ ; 

For  an  arbitrary  open  set  w  £  IR™  ,  m  =  1,2,  we  use  the  terminology 

H  (id)  to  denote  the  set  of  functions  that  together  with  all  their  derivatives 

of  order  <  k  are  square  integrable  in  u  .  Spaces  with  noninteger  indices 

are  defined  by  complex  interpolation.  If  g  is  an  analytic  curve  with 

(g(s)  :  s  €  (0,1)}  =  T  ,  then  HC(r)  denotes  those  functions  u  for  which 

u°g  5  Ht((0,l))  .  The  norm  on  HC  (u>)  is  denoted  j*  and  similarly  for 

t , 

t 


the  space  H  (2)  . 

Let  ft  be  an  analytic  curvilinear  polygon.  Ue  recall  that  _\g  * 
is  a  corner  if  near  it  ft  is  analytically  isomorphic  to  either  {jt  :  x^>0  and 
x?>0}  or  (x  :  x^>0  or  x2>(^^  •  Let  x(s)  »  s  *  [  —  5,6]  ,  5>0  ,  be  a 

parametrization  of  3ft  near  Xq  =  x,(0)  according  to  arclength  (i» e.  the  arclength 
from  x(s)  to  is  ]s|  ,  and  x(s)  and  x(-s)  lie  on  opposite  sides  of  x^). 

If  <p  is  a  function  on  3ft  ,  we  set 

r5,  ,2  -1 
I  (;)  =  j  j$(x(s))-d(x(-s))|  s  ‘'ds  . 

*0  J  0 

Let  x.  ,  1  <  i  <  M  ,  be  all  the  corners  of  32  ,  and  ]et  2.  =  g.(s)  :  s  r  (0. 

— 1  11 

where  g^  ,  1  $  i  <  M  ,  are  the  analytic  boundary  curves  connecting  these  corners. 


We  define 


H1/2(3ft)  =  {<J>  :  $j  €  H1/2(r.)  ,  1  <  i  <  M  and 

i  1 


viw.i  norm 


I  (<J>)  <  00 
x . 

— l 


M 


1  <  i  S  '[} 


’!*'!-l/2  “(I  ["*«'!;  +!(•)]) 

Hi/Z  i-i  |,r.  x- 


1/2 


— x 


1  ~l/2 

There  is  a  continuous,  surjective  trace  operator  from  H  (ft)  to  H  (3ft)  [!]• 


’  ■  '.•*»*.***  «.'•  .  * 


■.••'.s'. .  .  •  ^ 


*  •  * 


12 


the  statement  "u  =  <j>  on  3fl"  should  always  be  interpreted  in  the  sense  of  this  trace. 

00 

For  any  y  €  L  with  0  <  Yq  -  Y (x)  we  denote  by  the  correspond¬ 

ing  operator 


L^u  =  7*(yVu) 


The  equation 


L  u  =  0  in  H  ,  u  =  $  on  3* 


1  -]/o 

has  a  (unique)  solution  u  €  H  (Q)  exactly  if  €  H  (3fi)  ;  the  energy 

/j.  y|7uj2dx  is  a  continuous  quadratic  form  on  H^^(3K)  .  Whenever  u  £  K^"(.2) 
2 

and  L  u  €  L  (Q)  ,  the  conormal  derivative 


Y  ~  ^  [H1/2OSl)]*  (the  dual  of  H1/2) 
a  V 


is  defined  by 


f  3u  ,  f  ( 

y  y  ds  =  y Vu  •  ■  v  d>c  +  j  L  u'V  dx  , 

•'an  if!  i  a  Y 


where  v  €  H  (n)  is  any  function  satisfying 


v  =  on  ,  l|v;|  n  <  C^'l  .  n  . 

’  H  ^ 

The  Dirichlet-to-Neumann  operator  A  :  u!.^  -»•  y  maps  H^“(3ft)  to 

y  '  du  dv  cfi: 

- 1/2  * 

[H  “(3C!)]  ,  and  its  image  is  the  subspace  annihilated  by  the  constant  functions. 

We  turn  next  to  the  Runge  Approximation  Property.  Let  m  be  a  subdomain 
of  Q  and  consider  a  solution  of  L^u  *  0  in  j  .  It  is  generally  not 
possible  to  extend  u  to  a  solution  of  L  u  *  0  in  the  full  domain  P.  ;  that 
would  correspond  to  solving  the  Cauchy  problem  for  this  elliptic  equation.  It 


a;. 


is  possible,  however,  to  find  an  approximate  extension"  in  a  certain  sense. 
This  follows  from  results  of  Lax  [9]  and  Malgrange  Lll];  we  give  a  complete 
proof  for  the  reader's  convenience. 

Lemma  2.  (The  Runge  Approximation  Property)  Let  u  be  a  C°°  domain 
contained  _in  the  analytic  curvilinear  polygon  ?.  . 

and  such  that  each  connected  component  of  51\a)  has  a_  boundary  curve  in  common 
with  351  .  Let  0  <  y^  <  y(x)  b£  piecewise  analytic  on  55  and  assume  that 
u  €  H^w)  satisfies 


L  u  *  0  in  u  . 

Y 

Given  any  compact  subset  K  c  ui  and  any  e  >  0  there  exists  U  i  H^(51)  such 


L  U  *  0  in  *1  ,  and 
Y 


f  2 

I  |v(u-u) !  dx  <  e 
>K 


2 

Proof:  It  will  suffice  to  show  that  f  |U-u|^dx  can  be  made  as  small  as 

-  - - —  *  it\  — 


desired,  since 


1 7  (U-u)  |2dx  <  C  |lT-u|2dx 
•'K  ■'w 


as  a  consequence  of  the  identity  L^(U-u)  “0  in  u>  .  So  we  must  show  that 


H  «  { u  :  u-U|  ,U€H1(52),L  U»0  in  51} 


is  dense  in 


(u  :  u€L  (oj)  ,  L^u  “0  in  u>} 


vv-v-'.-; :v- •  -y 

C  v"  C  »■*  •-V'*-  C  v  i-  ■-  --  --  •-  •-  •-  '■  V  A  •-  •v  •-  v  ^  -V  •.  .  V _ _  -  •-  w  S-  «-• 


in  the  norm  of  L  (id)  . 


If  not  then  there  must  be  some  u*  €  L  (w)\{0)  satisfying 


L  u  =  0  in  w 
Y 


(7) 


j  u*v  dx  =  0 


for  all  v  €  H 


We  shall  derive  a  contradiction  from  this.  Let  $  €  H  (ft)  solve 


(8) 


L  $  »  U  in  ft  ,  $  -  0  on  3ft  , 

Y 


where  U*  »  u*  in  w  ,  U*  *  0  in  ft\tu  .  By  (7) 


0  -  u*V  dx 


L  dx 

Y  ~ 


I  Y  T7T‘V  ds 


3ft 


3v 


for  any  V  €  H^(ft)  with  L  V  =  0  in  ft  ;  therefore  y  *  0  on  3ft  . 

Y  3  v 


We  recall  from  (8)  that  $  also  vanishes  at  3ft  .  Since  y  is  piecewise  analytic 
★ 


on  ft  ,  U  *  0  in  ft\w  ,  and  each  component  of  ft\u>  has  a  boundary  curve  in 
common  with  3ft  ,  it  follows  by  repeated  application  of  Holmgren's  Uniqueness 
Theorem  (see  e.g.  [A])  that 


$  =  0  in  ft\m  . 


3$ 


Therefore  =  y  —  =0  on  3w  ,  and  as  a  consequence 

o  V 


r  *  2  r 

|  |u  |  dx  =  I  L  ^’u^djc  =  0 
J 7,.,  Y 


This  contradicts  the  hypothesis  that  u  was  nonzero,  and  the  proof  is  complete. 


Remark  2 :  It  is  not  essential  that  w  be  smoothly  bounded  or  that  ft  be  an 
analytic  curvilinear  polygon.  The  hypothesis  that  y  be  piecewise  analytic  can 
also  be  relaxed,  by  using  Aronszajn's  unique  continuation  theorem  instead  of 
Holmgren's.  The  version  stated  here,  however,  is  exactly  what  will  be  needed 
in  the  proof  of  Lemma  4. 


16 


4.  Energy  estimates 

The  following  is  a  slight  reformulation  of  a  result  in  [6]. 

OD 

Lemma  3.  Let  u>  be  a  bounded  C  domain,  Xq  a  point  on  3oj  ,  and  D  a 

_  30 

neighborhood  of  x^  relative  to  w  .  Assume  that  y  (  L  (w)  ,  0  <  J  y (x)  , 

00  1 
and  that  y  f  C  (D)  .  Given  1^0  and  e  >  0  there  exists  u  €  H  (w) 

such  that 

(i)  Lu=0  in  o 

Y 

(ii)  I  )Vu|^dx.  <  e  I  p^|Vu|^dx 

w\D  'D 

with  o  (x)  =  dist  (jc,  3u)  . 

Proof :  For  any  fixed  integer  M  >,  0  ,  Lemma  1  of  [6]  constructs  a  sequence 

{$„}!!  ,  c  C  (3u)  such  that 
N  N=1  — 

(9a)  !!*«:!,  «  C  NC  ,  t  5  -  M  , 

7+t,3w 

(9b)  UH!I  i  -  1 

J, 

(9c)  supp  4-  {Xq}  as  N  -*■  *  . 

We  review  the  construction  when  Xq  is  the  origin  and  *  {_x  :  x.,  =  0}  (to 

00 

which  the  general  case  is  easily  reduced).  Let  be  a  C  function  on  1R  , 

not  identically  zero,  with 

supp  4>  c  [-1,1  ]  , 

+1  k 

s  ^(s)ds  =  0  for  integers  k,  0  <  k  <  M  -  1  . 

-1 


4S  ,V 


If  )  *  ^(Nx.)  ,  one  easily  verifies  that  cb 

N  1  1  y  N 


V  ! 1 ,  has 
‘  7.3“ 


properties  (9a-c). 


Now  consider  the  solution  ul,  of 


\UN  =  °  in  “  >  “n  *  *N 


on  . 


Lennna  2  of  [6]  shows  that 


5  C  N' 


while  by  Lemma  3  of  [6] 


pe|VuN|2dx  5  C6  N"(2+6^  ,  Cfi  >  0  , 


for  any  5  >  0  .  A  combination  of  these  gives 


l_  1 2  _  „(2+6)«-2M  \  l,  ,2, 

lVuNl  d2S  *  c5  N  j  P  |VuN|  dx 


By  choosing  M  >  i  and  0  <  5  <  -j - 2  we  obtain  (2+6)  £  *  2M  <  0  t  and  hence 

V 


(2+SK-2M  .  c 

6 


for  sufficiently  large  N  .  The  corresponding  satisfies  (i)  by 

definition  and  satisfies  (ii)  as  a  consequence  of  (10)  and  (11),  using  the 
fact  that  J  p^|VuN|2dx  i  0  . 


IS 


The  next  lemma  combines  the  result  in  Lemma  3  with  the  Runge  Approxima¬ 
tion  Property  to  obtain  a  similar  energy  estimate  for  subdomains  and  with 
functions  that  are  "y-harmonic"  in  the  larger  domain.  For  our  later  app 1 i ca¬ 
tion  of  this  result  it  is  essential  that  the  subdomain  be  only  piecewise 

"smooth".  If  a)  is  an  analytic  curvilinear  polygon  then  we  shall  call 

3Cq  €  a  regular  boundary  point  (relative  to  y  )  if  there  exists  a  neighbor- 

2 

hood  B  of  x^  in  H  and  a  map  4>  such  that 

(i)  $  maps  B  injectively  onto  a  neighborhood 

2 

of  the  origin  in  , 

(ii)  both  maps  i>  and  0  ^  have  analytic  coordinate 
functions, 

(iii)  4>(B0u)  is  given  by  (x  :  x2>0}  D  $(B)  and 
y  is  analytic  on  Bflw  . 


Lemma  4.  Let 


be  two  analytic  curvilinear  polygons  with  ui  £  f< 


Let  hy  £  piecewise  ana  1  vr  ic  fun_c  t  i  on  on  with  0  q*  »(x1  .  Assume  that  Xq  €  3w 

is  regular  relat ive  to  y  and  that  x^  lies  in  the  unbounded  component  of 

2  _ 

1R'’  \  w  .  Let  D  be_  a_  neighborhood  of  x^  relative  to  oj  .  Given  i  >.  0  and 

e  >  0  there  exists  U  €  H^(ft)  such  that 


(i)  L  U  =  0  in  ft 

y 


(ii)  :  |VU£dx  <  £  J  p^jTl'I^dx  , 

j  w\  D  '  D 


where  p(x)  =  dist(x,3w)  . 


•*.  ■* 


*.  s  *  .*•  .*•  .*•  »'»  s  .  *  I 


9 


Proof : 

We  may  without  loss  of  generality  assume  that  <U  c  ft  and 
simply  connected;  if  not  we  merely  enlarge  fi  and  extend  y  by 

CO 

piece  of  ft  .  Let  u'  c  !)  be  a  C  domain  so  that 

(i)  b )'  is  simply  connected  and  contains  w  , 

(ii)  the  boundary  of  oj *  coincides  with  the 
boundary  of  u  in  a  neighborhood  of  Xq  » 

(iii)  w\D  c  to'  . 

(See  Fig.  2.) 


Fig.  2 


Lemma 


Choosing  a  neighborhood  D'  c  D  of 
3  to  construct  u  €  H^(uj')  satisfying 


relative  to  oj '  , 


that  ft  is 
1  on  this  new 


we  apply 


(12) 


|  |Vu|  d x  <  6  (P ' ) y I Vu I  dx 

Jtu’VD'  'n» 


where  p'(x)  =  dist(x,  3u)')  and  5  >  0  is  arbitrary.  If  D'  is  small 
enough  we  have 


A 


P '  (x)  =  dist.  (x,  3ui' )  =  dist(x,  3cu)  =  p  (x) 


for  all  x  €  D'  ,  and  then  (12)  implies 


(13) 


JW\D 


Vu  I  dx  <  6 


p  ^ I Vu I ^dx 


1 D 


Since  the  inequality  in  (13)  is  strict,  there  is  a  compact  subset  c  interior (0) 
for  which 


(14) 


j  Vu (  dx  <  6 


u)\D 


p^  [  Vu [^dx  . 


Let  K  =  J\D  U  Kp  .  Since  it  is  a  compact  subset  of  u'  ,  Lemma  2  applies  to 
give  a  function  U  €  h\;T2)  satisfying  L^U  =  0  in  $  and 


J 


'  |v(U-u)|^dx  <  S  j  p^jvuj^dx  . 

^ADUKp  J* 


Assuming  for  simplicity  that  p  <  1  in  ,  we  conclude  that 


| V (U-u) j “dx  + 


u>\  D 


p*' j7(U-u)  1  “dx  <  5  [  o^lTUi^dx 


*4> 


X  2  /  * 

It  follows  using  the  inequality  -jjaj  i  |a-b|“  +  jb|“  that 


(T  -  5) 


p^jVuj^dJC  <  j  p^lvUJ“dx 


(13) 


(16) 


X 

2 


|vu|  dx  <  5 


p  |  Vu !  d_x  + 


|  Vu  I  d}{  . 


Inequalities  (14),  (15)  and  (16)  combine  to  give 


J^\D 


(VUpdx  <  46 


f  p  7 

j  p  j7u|  dx 

% 


86 

1-26 


p^|7U|2dx 


$ 


86 

1-26 


P^'|VU|2dx  . 
D 


We  get  the  desired  inequality  by  choosing  6  so  that  85/ (1-25) 


Let  I?  be  an  analytic  curvilinear  polygon.  Our  goal  is  to  show  that 


two  piecewise  analytic  conductivities  and  y?  that  produce  the  same 

boundary  measurements  are  necessarily  equal.  It  is  assumed  that  y  ,  y^  <z  !°°(n) 
with  0  <  Yq  <  y^ (x)  (i=l,2)  .  "Producing  the  same  boundary  measurements"  means 

that  the  maps 


-1/2 

Hx/~(30) 


[ h1^2 ( ar;)  ]* 


defined  by 


3u . 

A  =  Y .  -r —  where  L  u.  =  0  , 

Y .  x  3v  y .  l  ’ 

l  l 


$ 


are  the  same  for  i  =  1,2  .  Equivalently,  it  means  that  the  quadratic  forms 


Q  (♦) 
i 


V  !?U  I2dx 
in  1  1  “ 


are  the  same. 

Theorem  1.  Let  y^  ,  i  =  1,2  ,  b£  piecewise  analytic  functions  on  ft  with 
a  positive  lower  bound .  If 

Q  (0)  =  Q  (i) 

Y1  Y2 

.  -1/2 

tor  all  :  -  H  (32),  then 


Y1  = 


Proof :  K’e  shall  assume  that  Y ^  and  .  y0  are  piecewise  analytic  relative 

N 

to  the  same  covering  tw.)  ;  by  Lemma  1  this  represents  no  loss  of  generality 


,  and  we  order  the  cover  so 


Seeking  a  contradiction,  we  suppose  that  y^  t  Y 2 

that  {w.}.  ,  are  exactly  the  elements  where  y,  #  y„  . 

1  J-1  1  2 

For  any  analytic  curvilinear  polygon  12'  ,  the  outer  boundary  is  that 
part  of  312'  which  lies  in  the  unbounded  component  of  .  We  claim  that 

there  is  an  analytic  curvilinear  polygon  12*  c  12  with  the  following  properties: 


(17a) 


12'  contains  all  u  ,  1  £  j  £  K 


(17b) 


The  outer  boundary  of  12'  has  an  analytic 
curve  in  common  with  3u>  for  some  j  , 

1  <  j  <  K  .. 


Indeed,  if  3oj^  has  a  curve  in  common  with  the  outer  boundary  of  12  for 

some  j  ,  1  <  j  $  K  ,  then  it  suffices  to  take  12'  =  12  .  Otherwise  we  choose 
K  _ 

x'  €  3(  U  w.)  and  x"  on  the  outer  boundary  of  12  such  that  x.'  anc* 


j-1 


J 


are  connected  by  a  piecewise  analvtic  (e.g.  piecewise  linear)  curve  lying 

K  _ 

entirely  within  12x11  to.  ,  except  for  the  endpoints.  Wo  may  also  assume  that 


K 


j-1 


K 


3  (  U  to.)  is  an  analytic  curve  near  x'  >  anc*  that  :J  w.  lies  locally  on 
j-1  J  K  j-1  2 

one  side  of  3(  'J  io.)  .  By  excising  from  12  a  tube  along  the  curve  connecting 

j  =  l  j 

_x ’  and  x”  (with,  say,  piecewise  linear  boundaries)  we  obtain  an  analytic 
curvilinear  polygon  12'  with  the  desired  properties.  (See  Fig.  3). 


outer  boundary  of  S 


Shaded  area  marks 


For  such  ft' 


,  let  Wj  ,  1  ^  j  £  K  ,  be  a  piecewise  analytic  domain  sharing 
an  analytic  curve  F  with  the  outer  boundary  of  ft'  .  Since  y^  and  y^  are 
different  on  uk  ,  their  Taylor  expansions  must  differ  along  F  .  Therefore 
there  is  a  point  in  the  interior  of  F  and  an  integer  l  0  such  that 


|Yj_(x)-Y9(x)  |  i  c  [dist  (x,  3Q * )  ]  ^  ,  c  >  0  , 


in  an  ft1  -  neighborhood  D  of  Xq  •  Note  that  the  boundary  point  Xq  is 
regular  relative  to  y^  ,  i  =  1,2,  and  lies  in  the  unbounded  component  of 


a. 

1R  \f.’  . 


Writing  p  (x)  =  dist(x, 3ft')  ,  and  relabeling  if  necessary,  we  have  that 


(18) 


Y^ (x)  -  Y2(x)  5  C  p  (x)  in  D 


Lemma  4  with  w  =  ft'  and  y  =  y  yields  a  function  U  €  H  (ft)  such  that 


(19a) 


L  U  =  0  in  ft 
Y1 


(19b) 


ft'\  D 


jvU|2dx  <  e  [  p^(vu|2dx  , 


where  e  >  0  is  a  small  parameter  to  be  chosen  later.  Setting 


i‘  =  U|3fi  €  H1/2(3ft)  , 


we  have 


(20) 


Q  M 

Yl 


Y1|VU|  dx 


=  [[  +1  +[  Iy. |VU|2dx  . 

'■■'FA ft'  Jft*\D  1 D'  1 


Since  y^  =  y^  outside  ft  the  first  integral  equals 


1  Y0|VU|2dx  ; 


i  i  9 

JQ\Q’ 


the  third  integral  can  be  bounded  from  below  using  (18), 


Y,  [  VU  |  ^dx  >  j  y  |  VU  |  ^dx  +  c  I  p  '  |  7U  j  2dJC 
1  JD  2  '  D 


Substitution  into  (20)  yields,  after  deletion  of  the  second  integral 


r  2  (  -> 

Qv  W  5  y0|71'|  dx  +  Y->  1  VU  j  “dx  + 

Yi  Jn\n'  2  h  2 


c  p^|VU|2dx 


Y?|VU!2djc+  c  p^|VU|2dx  -  Y2  I  VU  | 

1  2  JD  JQ’\D 


Applying  (19b), 


Y o | VU |  dx  $  (sup  Yo )  | VU |  dx 

)  “  Jp.’\D 


&  2 

<  t  •  sup  y  p  |  VU  [  doc 

2 


If  e  <  c/sup  y 2  (c>0  is  the  constant  in  (21))  then  it  follows  that 


2 

Q  M  >  y9|vu|  die 
Yi  Jn  2 


Q  (v)  =  min  y7I'7v|  dx  , 
Y2  v€H1(H)-,n 


6.  Lavered  structure 


In  the  class  of  L  (or  even  C  )  coefficients  y  ,  ident if iability  by 
means  of  boundary  measurements  remains  an  open  problem.  It  is  natural  to 
expect,  however,  that  the  layered  case  -which  is  essentially  one  dimensional  - 
should  be  easier.  In  fact,  in  this  case  the  Dirichlet-  to  N’eumann-data  map 
determines  spectral  data  for  a  certain  potential,  by  means  of  analytic  continua¬ 
tion  in  Fourier  space.  One  dimensional  inverse  spectral  theory  is  rather  well 
understood,  and  it  will  allow  us  to  recover  a  conductivity  y  which  is  merely 
three  times  differentiable.  We  wish  to  chank.  D.  Stickler  for  suggesting  the 
use  of  analytic  continuation  to  pass  from  boundary  measurements  to  spectral  data. 

Throughout  this  section  2  will  be  the  infinite  strip  (x  :  -««x^<"  , 

2  ® 

0*x7<l;  £  K  and  y  -  y(x2)  will  be  in  L  ((0,1))  with  a  positive  lower 

l/'1  2  1 

bound.  For  any  pair  €  [H  OR)]  there  is  a  unique  u  €  H  (0)  such  that 

(23)  V*(y(x2)Vu)  *  0  in  Q  ,  u|x  ,  u|^  • 


The  associated  Dirichlet-  to  Neumann-data  map  is 


a  («,;) 

Y 


5u  |  3u  | 

Y  3x2  lx?a‘0  ’  Y  ^x2  x2“l 


-1/2  2 
6  [H  OR)  \ 


We  begin  by  observing  that 
2 

when  y  is  merely  C “  . 


A  determines 

Y 


y  and 


at  3D  even 


Lemma  5 :  If  y  and  y2  are  C  on 
and  first  derivat ives  agree  at  3D  . 


and  A  =  A  then  their  values 
Y1  Y2 


Proof:  We  shall  check  that  the  argument  of  [6]  applies.  First  of  all,  the 

2 

regularity  estimate  (7)  of  that  paper  holds  for  k  =  2,3  when  Y  is  C 
Therefore  formula  (14)  remains  valid  for  M  •  2  ,  and  the  proof  of  Lemma  3  works 


with  t  «  3  in  (23).  If  l  <  1  in  (26)  then  the  argument  on  pp.  296-7  applies 
using  M  *  2  ,  and  this  yields  the  desired  conclusion. 


Remark  3:  This  proof  does  not  use  the  layered  structure,  but  it  does  use  that 

2  -  IR2  .  For  (2  c  ]Rn  the  same  argument  requires  y  €  Cr  ,  r  >  y  .  Taking 

into  account  the  layered  structure  it  is  possible  to  reach  the  conclusion 

that  r .  =  '•  from  only  knowing  that  A  ($,i|>)  =  ..  (4>,^)  for  two  sets 

Y1  Y2 

of  Dirichlet-data  (<p , iJj)  (see  Remark  5). 


Our  next  task  is  to  separate  variables.  We  shall  write  v^  *  v^(£,x2) 
for  the  Fourier  transform  of  a  function  v(x^,x^)  with  "respect  to  its  first 
argument.  (If  v  only  depends  on  x^  then  we  write  v(£)  ;  no  transforms  will 
ever  be  taken  with  respect  to  . )  From  (23),  the  solution  of  L^u  =  0  with 
Dirichlet-data  (<J>,4>)  satisfies 


(24)  dx_^Y(x2)  dx~  ^ (^ >x2>  >  ~  C2y(x2)u1(C,x2)  =  0  ,  0  <  x2  <  1 

^U.O)  -  i(c)  ,  u1(C,l)  »  *(5) 


for  a.e.  C  €  R  . 

The  ODE  (24)  has  its  own  boundary  data  map  A  for  each  fixed  £  €  K 

Y  y  Si 

2  2 

It  takes  1R  to  1R“  ,  mapping  (a, 8)  to 


A  (a, 8)  *  (-yv  (0),yv  (1))  , 
Y » s 


where  v  €  H  ((0,1))  solves 


(yv')’  -  5  Yv  ■  0  in  (0,1)  ,  v(0)  =>  a  ,  v(l)  -  3 


The  two-dimensional  boundary  data  map  A  determines  this  A  r  for 

Y  Y » •» 

every  £  : 

OD 

Lemma  6:  If  y.  and  y„  are  L  with  a  positive  lower  bound  and  A  *  A 
12 

then  A  „  *  A  _  for  each  £  €  IR  . 

-  Y^S  Y2»C - 

Proof :  Choose  Dirichlet-data  0  and  0  such  that  $  and  0  ,  when  restricted 
to  some  interval  (a,b)  ,  vanish  at  most  on  a  set  of  measure  zero.  From  the 
hypothesis  that 


A  (0,0)  *  A  (0,0)  and  A  (0,0)  =  A  (0,0) 


it  follows  easily  that 


A  *  A  for  a.e.  £  €  (a,b) 

YX.C  Y2.^ 


Since  A  is  real-analytic  as  a  function  of  £  ,  the  same  identity  must 

Y  9  S* 


hold  for  all  £  *  R 


□ 


In  order  to  apply  inverse  spectral  theory,  we  must  transform  the  ODE 

(26)  to  a  more  convenient  form.  If  v  solves  (26)  then  w  =  y^v  solves 

(27)  -w"+qw+£^w  =0  in  (0,1)  ,  w(0)  =  y^(0)a  ,  w(l)  =  yii(l)6  i 


with 


(28) 


-  t(y’)2/y2  +  h”/y 


We  note  that 


w  '  =  y**v '  +  4y  . 


The  Dirichlet-to-Neumann  map  for  (27)  is 


A  :  (w(0) ,w(l) )  -  (-w'(O).w'(l))  ; 
H  >  s 


clearly  it  is  determined  by  A  _  and  the  values  of  y  and  y ’  at  0 

Y  *  s> 

and  1 . 

Inverse  spectral  theory  addresses  the  problem  of  determining  the 

potential  q(x)  given  knowledge  about  the  eigenfunctions  and  eigenvalues  of 
2  2 

-d  /dx  +  q  .  The  result  we  shall  use  is  this  one,  proved  in  (12]: 

Lemma  7 .  Given  two  potentials  q  €  ( [ 0 , 1  ] )  ,  i  =  1,2,  consider  the  operators 


dx 


with  homogeneous  Neumann  boundary  conditions.  Let  A^  and  be  the 

associated  eigenvalues  and  eigenfunctions,  with  the  normalization  4/^  (0)  =1  . 

If  A^  =  A^  and  4>^(1)  =  $^(1)  for  each  m  then  q,  =  q,  . 

—  m  ra  -  m  m  -  -  -  1  z 

It  remains  to  relate  the  boundary  data  maps  A^  ^  and  the  spectral  data  for 

q  .  To  this  end  we  consider  yet  another  map  M  _  ,  taking  Cauchv  data  at  0 

q.S 

to  Cauchy  data  at  1  for  the  ODE  (27).  Explicitly, 


M  (a, 6)  =  (w(l),w'(I))  where 

q » s. 


-w"+qw+^  w=  0  in  (0,1)  ,  w(0)  =  a  ,  w' (0)  =  5 


It  is  easy  to  check  that 


M 


.  -  M 

qL.c  q2>^ 


if  and  onlv  if  A  _  =  A 
ql 


qn » c  q2>£ 


for  any  q_^  of  the  form  (28)  and  any  £  €  1R  (that  q^  has  the  form  (28)  is 


Mit  wav  to  guarantor  that 


is  vi*  1  i  —  df  t  i nod  > .  For  the  purpose  of  identifying  y  th 


use  of  M  is  thus  equivalent  to  that  of  A  .  An  important  difference 

qO  q,£  K 

is  that  M  F  is  well-defined  for  complex  arguments  £  (and  arbitrary  q  €  L  ); 
S  »  1 

indeed  it  extends  to  a  holomorphic  function  in  the  entire  complex  plane.  The 
values  of  M  .  for  complex  £  are  therefore  determined  by  those  for  real 

M  9 

£  by  analytic  continuation.  It  is  easv  to  show  that  M  (with  complex  £  ) 

q.s 

determines  the  spectral  data  of  q  : 

Lemma  8:  Let  q^.  ,  i=l,2,  be  two  potentials  in  C1([0.1]>  .  If  M  _  =  M 

1  —  q1,£  q2,<» 

for  all  complex  £  ,  then  q^  =  * 


Proof :  Assume  that  X  ^  ^  and  4/^  is  an  eigenvalue  and  the  corresponding 

2  2 

eigenvector  for  the  operator  -d  /dx  +  q^  with  Neumann  boundary  conditions, 

-  U(1)]"  +  q1«(1)  -  X(1),(1>  =0  in  (0,1) 


4> (1)  ’  (0)  =  b(1)’(D  =  0  . 


As  in  Lemma  7,  we  may  take  the  normalization 


(0)  =  1 


2  (1) 

Choosing  n  £  £  so  that  n  =  -X  ,  and  noting  that 


«  „  u.°>  =  „  d,o)  , 


we  see  that  the  solution  of 


-w"fq  w-X '  yw  =  0  ,  w  (0)  =  1  ,  w' (0)  =  0 


additionally  satisfies 
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w(l)  -  '(1)  ,  w’(l)  -  0  . 

(1)  2  2 

Therefore  A  is  also  an  eigenvalue  of  -d~/dx“  +  q^  ,  w  tfie  corresP°nd" 
ing  eigenvector,  and  its  Cauchy  data  agrees  with  that  of  at  both  end¬ 

points.  A  similar  result  holds  by  symmetry  for  any  eigenvalue  and  the 
corresponding  eigenvector  for  q?  .  From  Lemma  7  we  now  conclude  that  =  q? 


We  obtain  the  identif iability  of  a  C  ,  layered  conductivity  by  assembl¬ 
ing  these  results. 


Theorem  2:  Let  P.  be  the  strip  {x  :  -°°<x^<®,  0<x^<l}  .  If_  y^(x^)  and 

3 

Y,,(x„)  are  two  layered  conductivities  of  class  C  ([0,1])  ,  and  if  A  =  A 
“  “  Y1 
(or  equivalently  if_  Q  -  Q  )  ,  then  y  =  y0  . 

Y1  Y2  1 

Proof:  Let 


qi  *  ~  i(Yl)2/Yi +  k/Yi  •  1 "  1»2 


Lemmas  5  and  6  and  the  discussion  immediately  after  show  that  A  *  A  _ 

qj.5  q2*s 

for  all  £  €  ]R  ,  and  we  have  explained  why  this  implies  M  “  M  for 

>  L  q2 »  l 

all  £  €  <E  .  Therefore  q^  ®  q^  as  a  consequence  of  Lemma  8.  For  fixed  q^ 
(29)  is  a  second-order  differential  equation  for  y^  .  Both  y's  have  the 
same  Cauchy  data,  using  Lemma  5  again,  so  they  must  be  equal. 


Remark  4 :  One  could  use  other  inverse  spectral  theorems  in  place  of  Lemma  7, 

with  minor  modifications  of  the  argument.  The  one  proved  in  [5]  requires  less 

2 

regularity  for  q^  ,  and  it  leads  to  a  proof  of  Theorem  2  with  y^  €  C  ([0,1])  . 
We  believe  that  the  regularity  hypotheses  on  y  are  an  artifact  of  the  method 


.  *,  ■,  *»  %  *,  m‘.  •  t*.  - .  ’  • , " •  _  ~  .  *  .  •  .*• ,  • .  ■  .  ■ .  * .  ■  _» * .  -  ..  *_«  *  ■"*  *  *  _■  *_«  *_•  *..**.. *_* 


of  proof,  and  not  intrinsic  to  the  problem. 


□ 

Remark  5:  Wt-  nay  i  <>n*- 1  tide  that  Y.  =  i based  an  a:i  (apparent)  weaker 

-  i  i. 

assumption  than  '  =  Tt  suffices  to  know  that 

*  1  ?2 

A  (4>,0)  =  A  ($,0)  and 

1 1  y2 

A  (0,w)  =  A  (0,'i<) 

Y1  y2 

for  two  specific  pairs  of  boundary  data  ($,0)  and  ( 0 , tp )  satisfying  the 
condition: 

there  exists  an  interval  (a,b),  a  <  b  , 
such  that  representatives  of  $  and  , 
when  restricted  to  (a,b),  vanish  at  most 
on  a  set  of  measure  zero. 

Indeed,  with  this  assumption,  it  follows  directly  from  the  proof  of  Lemma  6 

that  A  =  A  for  each  £  S  ]R  .  An  energy  argument,  similar  in  spirit 

Y1  ’  ’  y2  ’  ^ 

to  that  in  [6]  (or  that  in  the  proof  of  Theorem  1)  but  applied  to  the  equations 

(YjV' ) '  -  C"y1v  -  0  , 

with  i  sufficiently  large,  now  shows  that  y. (x)  =  yn(\)  and  y! (x)  =  Y^(x) 

at  x  =  0,1  .  From  the  discussion  following  Lemma  6  it  follows  that  A  =  A 

ql’ •' 

and  the  rest  of  the  proof  of  the  fact  that  y,  =  y.  proceeds  exactly  as  before 
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7 .  Convergence  of  a  reconstruction  algorithm 

The  numerical  reconstruction  of  a  conductivity  y  involves  its  estima¬ 
tion  in  some  finite  dimensional  space,  using  partial  information  about  the 
boundary  measurements.  We  shall  formulate  a  simple  algorithm  of  this  type,  and 

prove  its  convergence  as  a  consequence  of  the  identif iabil ity  results  in  [6]. 

N 

If  the  available  data  are  the  energies  {Q  (?  )}  ,  corresponding  to 

Y  m  m=l 

some  known  boundary  values  •  ^2  ’  ‘ ’  an<^  ^  Y  is  to  be  approximated  in  a 

00 

finite  dimensional  subspace  of  L  (Q)  ,  then  a  natural  approach  would 

seem  to  be 

find  y„  €  V„  H { 7  :  0<  ;(x) }  such  that 

-  N  N  0  — - 

max  |Q  (i>  )-Q  ($  )!/;]$>  j|?  is"minimal"  . 

l*m«N  YN  ra  Y  m  mi,3fi - 


For  this  to  succeed,  it  is  of  course  necessary  that  y  be  approximated  well 
by  some  element  of  .  That  ip  easily  arranged,  but  it  is  not  enough:  since 

the  problem  is  ill-posed,  we  must  somehow  make  sure  that  there  is  a  minimal  value 
.j:.J  that  •,  can  not  deviate  wildly  from  .  A  natural  wav  to  achieve  this  is 

■  ’  choose  \  within  a  fixed  compact  set  K  (relative  to  a  topology  in  which  we 

expet t  conver gence ) .  Thus  modified,  our  method  is  to 


(30) 


find  y»,  (  V„  fl  K  such  that 
-  N  N  -  - 


2 

max  !q  (}  )-Q  (?  )|/ljb  ■; ,  is  minimal 
l<m<N  YN  m  Y  m  - 


Finding  a  K  which  ensures  convergence,  is  not  overly  conservative, 

and  has  a  simple  description  is  in  general  nontrivial.  The  method  (30)  could 

be  called  semidiscrete,  since  the  evaluation  of  0  (t  )  still  requires  the 
-  -y  m 


-V*  y 

.•.sc 


%  , 


.V  . 

-o  o 1 — 


l  - 
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solution  of  an  elliptic  boundary  value  problem.  The  compactness  of  K  ensures 

that  the  sequence  {y.,}.,  .  has  a  subsequence  converging  to  some  y*  Z  K  .  If 

N  N=1 

00 

the  spaces  V.  ,  the  set  K  ,  and  the  sequence  {$>  }  ,  are  chosen  properly 

N  m  ra—  i 

then  y*  will  have  the  same  boundary  measurements  as  y  .  If  y  is  identifiable 
within  K  by  means  of  boundary  measurements  then  y*  =  y  and  the  method  converges. 

We  give  a  concrete  example,  including  the  details  of  the  convergence 
proof  just  sketched,  for  y  which  may  be  extended  holomorphically  to  a 


(complex)  ball 


=  :  ]z-x-!<R) 


W  =  l^€4:  •'  ‘-"V 


about  some  x.  €  2  ,  with  2  c  Bn(x_)  .  In  other  words,  we  suppose  that  y  has 
— -0  K  u 


a  power  series 


(x)  ■  l  aa(x-x0)ot  ,  x  £  P. 


sup  R 1 J  i a  I 


(In  the  last  two  expressions  a  ranges  over 


:  W 


all  .'-tuples  of  nonnec.at  ive  integers.) 


The  •‘unctions  •  t  .  are  assumed  to  be  dense  in  H‘O.1)  (  ..  is  a  C 
m  m=  1 

d.  rail > ,  and  we  choose 


VN  =  {p  :  p  is  a  polynomial  of  degree  £  d^ : 


where  ^jj^N~1  is  a  sequence  converging  to  00  as  N  00  .  The  compact  set 
K  consists  of  those  analytic  functions  ;  on  ..  tor  which 


—  min  (x)  i 


;  ;  C  ma  x  •.  ( x ) 
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sup  -V  r'3' |Da;(v  ) |  £  C  sup  R 1 3 ' | a  | 

Ct  •  ~~\J  ct 

a  a 


for  some  fixed  C  >  1  . 

It  is  easy  to  see  that  y  can  be  approximated  well  in  f!  K  ,  Indeed, 


when  N  is  large  enough 


=  )  a  (x-x_)  '  VM  n  K 

N  |  V  a - 0  N 

fa  I <dN 


and  this  -*■  y  in  L  (ft)  as  N  •  We  know  that 


Q?(<l>)-Q  (♦)  I  $  Djlc-yil 


l  (ft)  k.dn 


whenever  ?  and  y  are  uniformly  bounded  away  from  0  and  00  .  Therefore 


min  max  [Q  (<fr  )-Q  (<t>  )  |  /  ji  i>  !|,  r  -*■  0 


l<?ra<N 


C  tn  y  m  m  diw 


as  N  -►  ®  ,  and  so  the  y^  satisfy 


max  IQy  ao 

lgnsS  YN  m  Y  ra  ra'5’3n 


as  N  -*■  *  . 


3y  the  definition  of  K  , 


YN  =  I  a3*  (x-Xq)* 


vi  tn  a 


0  for  { ct  |  >  d^,  and 


_  01 

sup  R 

a,N 


£  C  sup  R  a 
1  a 

ct 


Therefore  there  is  a  subsequence  for  which 


(v 

a  b  as  k  -*■  00  .  for  each  a 

a  a 


and 
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sup  R^  b  $  C  sup  R^  a 
a  a 

a  a 


It  follows  ‘that 


l  b  (x-x  )a 
L  a  0 


uniformly  in  every  ball 


<  x 


x-x,. 


R'  ) 


R'  ■"  R  ,  and  consequently 


(35) 


Ym  -*•  Y*  in  L°°(a)  . 

Nk 


From  (33)  and  (35)  we  get  that 

Q  *(4  )  =  Q  ( 4  )  m=l,2, . . .  , 

Y*  m  Y  m 

and  thus 


( 36) 


v =  0 


since 


{*  }  is  dense  in  H  (3i!) 
m  m=l 


The 


using  ( 3-i  )  ,  (3"'»  and  the  fact  that 


'  N, 


•m  1  36)  and  tin  ident  i  ihil  it  v  «•>  f  an  tlvt  in  con due  t  i  v  i  t  ies  ,  proved  in  [6] 


tr 


follows  that 


The  nrocedin,’ 


has  a  convergent  subsequence  u-ndme 


argument  shows  not  only  that  ^j^N-1 

to  y  hut  also  that  any  convergent  subsequence  of  {y  }°°  necessar ilv 

N  N=1 

tends  be  y  ,  This  implies  that  the  entire  sequence  converges  to  y 


Remark  6:  The  case  when  y  is  merely  analytic  on  jj  can  be  treated 
similarly.  Let  H  be  a  compact  set  such  that 


H  c  interior(H)  ,  and 


y  extends  analytically  to  a 
neighborhood  of  H  . 

The  latter  condition  assures  that 


for  some  >  0  ,  where  the  supremum  is  now  over  points  x  ^  H  and  integer 
two-tuples  a  5  0  .  If  (31)  is  replaced  by 

sup  RHCt '  I  D°^  I  5  c  SUP  A-  r!0,  ^  |  (x)  I 

in  the  definition  of  K  (the  supremum  again  being  over  x  €  H  and  a  >  0) 
then  the  argument  can  proceed  essentially  as  before. 

Remark  7 :  A  similar  semidiscrete  method  can  be  formulated  for  the  layered 
case.  By  Theorem  2,  convergence  can  be  assured  by  choosing  K  so  that  four 
derivatives  of  y  remain  uniformly  bounded  on  2.  . 

Remark  8:  For  piecewise  analytic  y  we  know  ident if iability  from  Theorem  1,  but  it 
is  not  so  clear  now  to  define  K  .  If  a  cover  relative  to  which  K  is  piecewise 
analytic  is  known,  then  one  can  proceed  as  before  by  controlling  y^T  separatelv 
on  each  piece.  If  the  cover  is  unknown,  however,  it  seems  difficult  to  find  a 
choice  of  K  which  assures  convergence. 
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